GEOMETRIC REGULARIZATION ON RIEMANNIAN AND 
LORENTZIAN MANIFOLDS 
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^.J Abstract. We investigate regularizations of distributional sections of vector 

J^ bundles by means of nets of smooth sections that preserve the main regular- 

ity properties of the original distributions (singular support, wavefront set, 
Sobolev regularity). The underlying regularization mechanism is based on 
^^v ■ functional calculus of elliptic operators writh finite speed of propagation with 



< 

-)— » 

a- 



o 



L - 



^: 



respect to a complete Riemannian metric. As an application we consider the 
interplay between the wave equation on a Lorentzian manifold and correspond- 
ing Riemannian regularizations, and under additional regularity assumptions 
we derive bounds on the rate of convergence of their commutator. We also 
show that the restriction to underlying space-like foliations behaves well with 
respect to these regularizations. 



1. Introduction 



^Sl ■ We consider regularization processes to smooth out distributions on Riemannian 

^ I and globally hyperbolic Lorentzian manifolds and investigate their compatibility 

^^ ■ with the wave-equation. The Riemannian setting has been addressed in [13' and 

\^ I we include an introduction to this approach. For a globally hyperbolic manifold we 

^—N ■ pick a splitting of the metric as obtained by [2] which provides us with a globally 

defined time function and a foliation by space-like hypersurfaces. The associated 
Riemannian metric naturally allows us to construct regularization processes on 
distributional sections of tensor-bundles and differential forms. We show that these 
regularizations interact nicely with the wave-equation on the Lorentzian manifold 
and with the foliation provided by the metric splitting. 

By a regularization process we mean a net of smoothing operators that assigns a 
^\f ' net of approximating smooth functions (or sections in a vector bundle) to any given 

distribution (or distributional section). We are interested in preserving a maximal 
set of regularity properties of the distribution in this process (support, singular 
support, wavefront set, Sobolev regularity). To assign such properties to approxi- 
mating nets of smooth objects we employ the language of algebras of generalized 
functions in the sense of Colombeau ([9l [TOl 131] EQ]). In this approach any such 
regularization process provides an embedding of distributions into a space G{M) of 
nonlinear generalized functions (given by a quotient construction on spaces of ap- 
proximating nets). This process preserves the regularity and singularity structure 
of the distributions as described in Section [2] below. The main interest in employing 
these regularizations is in studying non-smooth curved space-times. 
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Our approach is motivated by work on wave equations on non-smooth curved 
space-times. C.J.S. Clarke (|6l[7l[8]) suggested to study physical fields for under- 
standing the singularity structure of the space-time itself, i.e., to consider wave 
equations in low regularity. His work gave rise to a number of further studies. 
In particular, in the framework of generalized functions this line of research was 
pursued in [34j [191 133] . For a more detailed introduction we refer to ^25] . 

Thus let us consider a model situation and analyze the generalized solutions 
to wave equations corresponding to a generalized Lorentz metric g on the smooth 
manifold M (cf. [IDl HZl HH])) i-e., a global Cauchy problem with initial data a, 6 G 
Ty'{S) on a suitable initial value surface S in the form 

DgU^O, u|5= i(a), Vnw|s= '-(fe), 

where i is assumed to be an embedding of distributions into the algebra of gener- 
alized functions. 

Suppose that the generalized Lorentzian structure induced by g allows a splitting 
of the wave operator in the form Dg = 7 • (Dg -|- Q), where g is a smooth "back- 
ground" metric, 7 is a positive measure or a strictly positive generalized function, 
and Q is a partial differential operator with distributions or generalized functions 
as coefhcients, but whose coefficient singularities are concentrated in certain space- 
time regions. A simple example is a generalized Robertson- Walker space-time with 
g = ~dP + {1 + fi{t))h on M ~ B. X S , where h is a smooth Riemannian metric on 
the smooth manifold S and fi is a nonnegative generalized function. In this case 
g = —dt^+h can serve as a smooth background metric and the above splitting would 
involve 7 = 1/(1 + fJ-) and Q = Q{t,dt), which acts only on the one-dimensional 
factor. 

If we impose the splitting assumption as above, then we have for any generalized 
function u G G{M) that the equation Dg u = is equivalent to Dg u + Qu = 0. Let 
u E G{M) be a solution to the Cauchy problem stated above. In attempting to 
extract its distributional aspects or assess its singularity structure we could adopt 
the following strategy. Let v be the distributional solution (assuming that it exists 
and is unique) to 

and put w :~ l{v) — w G G{M) as a comparison of the generalized function u 
with its "background distributional aspect" . Now, if the embedding l commutes 
with restriction to S and with Dg, then simple manipulations allow us to draw the 
conclusion that w satisfies the following Cauchy problem: 

DgW ^ Qu, w\s^O, Vnw|s = 0. 

Thus, w satisfies an inhomogeneous Cauchy problem corresponding to the smooth 
background metric g and with generalized functions occurring only on the right- 
hand side of the equation. In particular, w will be smooth in regions where null 
geodesies corresponding to g emanating from S do not intersect the C°°-singular 
support of Qu. 

In the more general situation of a globally hyperbolic Lorentz manifold provided 
with a foliation into spacelike hypersurfaces via a splitting of the metric we shall 
construct an embedding t and show that the embedding almost commutes with 
Dg and with restriction to space-like slices of the foliation. The extent of failure 
to commute is roughly measured by the divergence of the volume element on the 



slices along the time direction. Although we state our results in the scalar case, the 
extension to the case of wave-equations on differential forms is obvious. 

2. Regularization of distributions on complete Riemannian manifolds 

In this section we describe a geometric regularization procedure for distributions 
(or distributional sections of vector bundles) on complete Riemannian manifolds. 
This procedure encodes regularity and singularity features in terms of asymptotic 
behavior. Our approach is based on |13j . to which we refer for further details. 

We first fix some notations from the theory of distributions. Let M be an 
orientable complete Riemannian manifold of dimension n with Riemannian metric 
g. The space 'D'{M) of Schwartz distributions on M is the dual of the space 
n"(M) of compactly supported n-forms on M. Further, I?(Af) is the space of 
smooth compactly supported functions on M. We identify T>{M) with J7"(M) via 
f f-^ f ■ dg, with dg the Riemannian volume form induced by g. Thus I?'(M) can 
be viewed as the dual space of 'D{M). We consider L\^^{M) (hence in particular 
C°°(M)) a subspace of M via / i— > [(/p i— >■ J^^^ ffdg]. If -E is a vector bundle over 
M then I?'(M : E), the space of £'- valued distributions on M is defined by 'D'{M : 
E) = V'{M) (g)c=o(M) r°°(Af : E), with r°°(Af : E) the space of smooth sections 
of E. By £'{M : E) we denote the space of compactly supported distributional 
sections of E. We shall assume that E is endowed with a Hermitian inner product 
so that the distributional sections 'D'{M : E) oi E can be identified with the dual 
r^(M : E)' of the space of compactly supported smooth sections. Regularity 
properties of distributions we are interested in are, in particular, encoded in the 
singular support singsupp(w), the wavefront set WF(ii;), and the Sobolev regularity 
of any given w € 'D'{M). 

In order to be able to track these regularity properties in terms of regulariza- 
tions of a given distribution we need a conceptual framework that allows to assign 
geometrical and analytical properties to regularizations, i.e., to nets of smooth func- 
tions. Algebras of generalized functions in the sense of J. F. Colombeau ( UHl ISI] ) 
indeed provide a well-developed theory of this kind and below we shall briefly review 
some basic definitions, based mainly on [201IT6] . 

The basic idea of Colombeau's approach is to assign to any given locally convex 
space X a space Qx of generalized functions as follows. We define the space Mx 
of moderate nets to consist of those maps e >—>■ Xe (e € / := (0, 1]) that are smooth 
and for any seminorm p oi X satisfy p{Xg) — 0{e^) for some integer A^ as e — >■ 0. 
Similarly we call a net Xe negligible if for all seminorms p and all integers A^ the 
asymptotic relation p{Xg) ^ 0{e^) holds. The space Qx is defined as the quotient 
of the space of all moderate nets by the space Afx of all negligible nets. The 
class represented by a net x^ shall be denoted by [x^]. In case Af is a smooth 
manifold we call G{M) := Gc°°{m) the standard (special) Colombeau algebra of 
generalized functions on Af ([5J [T^l [2D])- If -E — >■ Af is a vector bundle then we set 
Q{M : E) :— Gr^{M:E)- For E = C the space Gc inherits a ring structure from C. 
It is therefore called the space of generalized numbers, denoted by C. Every space 
Ge is naturally a C-module, and hence is called the C-module associated with E 

(nsi). 

Similarly we shall also consider the subspace Gx of regular elements of Gx, 
defined as the space of all those elements of Gx that can be represented by a net 
Xg for which one can find an integer A^ in the above relations independent of any 



4 SHANTANU DAVE, GUNTHER HORMANN, AND MICHAEL KUNZINGER 

seminorm p on X . Thus Q^ is the space of uniformly controlled asymptotics. Again 
we set g°°{M) := G^^^^m)- I^ regularity theory, g°°{M) is the analogue oiC°°{M) 
in the theory of distributions ( |31[ [Ml IHl 113 )• This is based on the fundamental 
resuh ([31]) that for open subsets fl of R", L{V'{n)) D g°°{n) = C°°(17), where t is 
the standard embedding of V into G via convolution with a standard mollifier. We 
will introduce further notions of regularity theory based on G°° below. 

The assignment X i— > Gx is obviously functorial, so that any continuous linear 
map : X — > y naturally induces a map (/)* : Gx -^ Gy- Thus a smooth map 
f : M ^ N gives rise to a pull-back /* : G{M) -^ G{M). In particular G{M) 
defines a fine sheaf of algebras and similarly G{M : E) is a fine sheaf of G{M)- 
modules. 

We finally introduce the notion of association, which provides a concept of weak 
equality between elements of Colombeau spaces as well as a standardized way of 
assigning distributional limits to certain Colombeau generalized functions. For u, 
V G G{M : E) we say that u is associated (or weakly equal) tov.uK, v \iuir — Ve — >■ 
in 'D'{M : E) for (some, hence any) representatives of w, v. We say that u possesses 
w e 2?' (A/ : £■) as a distributional shadow if u^ — >• w; in I?'(M : E). 

To construct the desired regularization process on T)' {M : E) we need two pieces 
of data: 

(1) A Schwartz function F £ S^{M.) such that F = 1 near the origin. 

(2) An elliptic symmetric differential operator D on E such that the speed of 
propagation Cd, defined as 

(2.1) Cd = snp{\\aDix,m \ x G M, ||e|| = 1} 

is finite. Here gd denotes the principal symbol of D. Such an operator D 

will be referred to as admissible. 
As a consequence of the finite speed of propagation, the symmetric operator D is 
essentially self-adjoint ("231, 10.2.11]). Therefore the equation 

(2.2) —u^iDu m(.,0) = mo, 

has a unique solution for all times t for any initial datum uq G T^ (M : E) . Indeed 
it follows from functional calculus that e^^^uq is a solution, and uniqueness can be 
established using energy estimates. 

We note that functional calculus defines a map 

y{R) 3f^ f{D) e B{L\M : E)) 

where the operator ,f{D) can also be expressed by the Fourier inversion formula: 

(2.3) f{D) ■.= !-[ me^^'^ds 

(in the sense of strong operator convergence). Since the operator D is elliptic, it 
follows that /(D) has a smooth kernel by elliptic regularity. 

Let F G ^(M) with F = \ near the origin and set F^{x) := F{ex). Then we 
shall obtain the desired regularizing procedure from Fg{D). To do this we first take 
a closer look at singularity properties available in generalized functions. 

As has been mentioned above the space G{M : E) is a fine sheaf over M, hence 
provides a notion of support. Furthermore, the notion of wave-front set can be de- 
fined for generalized functions both locally and in invariant global terms analogous 



to distributions. For il C R" open we call u E G{^) tJ°°-microlocally regular at 
(xoj^o) G T*n \ if there exists some test function (p E 2?(rj) with (p{xo) = 1 and 
a conic neighborhood F C R" \ of ^o such that the Fourier transform F{ipu) is 
rapidly decreasing in F, i.e., there exists N such that for all Z, 

The generalized wave front set of u, WFg(u), is the complement of the set of points 
(xo, fo) where u is t/°°-microlocally regular. By [21], for any u E G{M), WFg(u) can 
naturally be viewed as a subset of T*M \ . It is equivalently defined as (cf. [ITj ) . 

WFg{u)= fl char(P) {Pe^%{M)). 

Here char(P) C T*M is the characteristic set of the order classical pseudodifferen- 
tial operator P. The singular support of u, singsuppg(M), is the complement of the 
maximal open set on which u E Q°° ■ It then follows that singsupp„(w) = pr]^(WFg). 
We are now ready to define a class of regularization procedures for distributional 
sections of a vector bundle. 

2.1. Regularizations. 

Definition 2.1. A parametrized family {T^)ir^i of properly supported smoothing 
operators (in the sense of [4, ch. 1.4]) is called an optimal regularization process if 

(1) The regularization of any compactly supported distributional section s E 
£'{M : E) is of moderate growth: For any continuous seminorm p on 
r°°(M : E), there exists some integer N such that 

piT.s) = O(e^) {e ^ 0). 

(2) The net (T^) is an approximate identity: for each s E <S"{M : E), 

limT^s^s inV'{M:E). 

e— >0 

(3) If u G T^{M : E) is a smooth compactly supported section of E then for 
all continuous seminorms p and given any integer m, 

p{T,u~u) = Ois"'). 

(4) The induced map lt ■ S'{M : E) -^ G{M : E) preserves support, and 
extends to a sheaf map 2?'(M : E) -^ Q{M : E) that satisfies, 

(2.4) lt{V'{M : E)) n g°°(M : E) = F°°(Af : E). 

This implies in particular that lt preserves singular support. 

(5) The map ut preserves wave-front sets in the sense that for any distribution 
s E V'{M : E) 

WF(u)=WFg(iT(u))- 

As mentioned already we shall use spectral properties of the elliptic differential 
operators to obtain regularizing processes satisfying the above conditions. 
Our main result in this section is the following. 

Theorem 2.2. Let F E ,y(R) he a Schwartz function such that F = \ near the 
origin. Let Fg{x) := F{ex). Given an admissible differential operator D, the family 
of operators {Ff;{D))i,^i provides an optimal regularization process in sense of Def 

^3 
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The essential idea of the proof boils down to the following two steps. 

First we prove Theorem 12.21 under the assumption that the underlying manifold 
M is compact. We then use finite propagation speed to extend the result to a 
general complete Riemannian manifold. 

Thus let us for the time being assume that M is a closed manifold. Then the 
space of smoothing operators ^^^(M : £') is a Frechet algebra. The functional 
calculus map defined by D induces a smooth map 4)d '■ ^(M) -^ ^^°(A/ : E), 
hence the estimate [T] holds because Fg, is a moderate net and F^{D) = (0£)),(Fe), 
where {4>d)* is the induced map on the asymptotic spaces ^,y(R) — > Q^-^[m)- 
Condition [2] then follows as the ring map 0c preserves approximate units. 

When M is compact the operator D has discrete spectrum and the spectrum of 
D^ satisfies Weyl's law, namely 

(2.5) Nd^{\) := #{A, G sp(I?2)| A, < A} - CA^^"^. 

The Weyl estimates in conjunction with the fact that F = 1 near the origin now 
provides the compatibility condition [3] (see dJ). 

To describe the regularity of a distributional section u E 'D'{M : E) we first note 
that any distribution provides a map between two Frechet spaces, namely between 
the smoothing operators ^^^(M : E) and the smooth sections T°"{M : E) by 
evaluation. More precisely, to any u S V'^M : E) we associate the map 

e„ : *^(M) ^ r°°(Af : ^) e„(r) := r(u). 

The mapping properties of the maps 0„ imply that if m ^ H^{M : E) for every 
k > t then given any J > 0, ||Fj(£))u||^2(j(_f.^) is not 0{e~^ +t+s^ j-ggg jjll Lemma 
7.4]). This direct description of Sobolev regularity is in fact stronger than condition 
4 and therefore, 

if^(i5)(I?'(A/ : E)) n g°°{M : E) = r°°(M : E)). 

The sheaf property of the embedding ip^ (£>) as well as the proof of the result for 
the more general case of a not necessarily compact complete manifold rely on the 
finite speed of propagation of D, to which we turn next. 

2.2. Finite speed of propagation. We shall reduce the computation of asymp- 
totics on a complete Riemannian manifold to certain compact manifolds obtained 
as doubles of suitable compact submanifolds with boundaries. If X is a compact 
manifold with boundary, a double of X, denoted here by DX is a closed manifold 
obtained by gluing two copies of X along the boundary dX (see, e.g., [26]). If X 
is a compact manifold with boundary embedded in a Riemannian manifold M of 
the same dimension and if U is an open subset of M such that U C interior(X), 
then there exists a Riemannian metric on DX such that the inclusion j : U ^^ DX 
becomes an isometric embedding. Moreover, for any vector bundle E -^ M there 
exists a corresponding vector bundle Ex — > DX such that Ex\u is canonically 
isomorphic to E\^. Also, there exists a symmetric elliptic operator Dx on Ex that 
coincides with D on U. 

Let u G £'{M : E) and fix a constant c > 0. Then by the Hopf Rinow theorem 
the open ball U := i32cCi3 (supp(m)) (with Cd as defined in (|2.1|) ) is relatively 
compact and therefore contained in a compact manifold with boundary X C A4 . 
By the above u can be identified with a distributional section of Ex -^ DX. 



Proposition 2.3. With assumptions on u,c and F as above F{D)u and F{Dx)u 
are both supported in U and 

F{D)u^F{Dx)u. 

Proof. The restrictions of D and Dx to tlie open set U coincide, hence uniqueness 
of solutions to (|2.2p implies that e*^^M and e'^^^u agree for s < c. Hence the claim 
follows from the Fourier Inversion Formula (12.31). D 



From these observations the proof of Th. 12.21 can be deduced along the following 
lines (cf. [El Sec. 4]): 

Let (j) G !){—€, c) satisfy = 1 in a neighborhood of and let u e £'{M : E). 
Then 

F,{D)u - ^y°° Hs)F,{s)e''''uds^j* [^ j^ (p{s)F,{sy'°''uds 

Based on this calculation we observe that the estimates required for Definition 
I2.1l in the general case of F^{D) follow from that of F^{Dx) which have already been 
established in case of the closed manifold DX. Furthermore the sheaf properties 
of i'Fe(D) can be obtained as a consequence of the calculations in Proposition 12.31 
The support of u coincides with the generalized support of [F^{D)u\. This implies 
that the embedding extends to a sheaf morphism ip^ : I?'(M : E) — > G{M : E). 
Preservation of wavefront sets is more involved. We refer to [13, Th. 3.10] for a 
complete proof of this property. 

Remark 2.4. In the scalar case, an alternative proof (not relying on the above 
doubling-technique) of Theorem 12.21 can be found in [T^, Section 3. 

2.3. Applications. As a first application consider E ~ /\* Af , the exterior bundle 
over a Riemannian manifold M and D = d -\- d* where d* is the Hodge adjoint 
of the de Rham differential d. Then D is symmetric, elliptic and has propagation 
speed Cd — 1 since a-£>{x,S,)^ ~ — ||^||^id, hence is essentially self-adjoint. Letting 
A :— D^ be the Laplace-Beltrami operator on fl*{M), it follows from functional 
calculus that cos(s-D) = cos(s-\/A) on L'^{M). 

Now let F be an even Schwartz function which has germ 1 at the origin. Then 

1 f°° 

F{^/A) = — / F{s) cos(sVA) ds 

as a Bochner integral in B{L'^{M)). 

Finally, let c > and pick an even test function (pc with support in (—2c, 2c) 
and such that 0c = 1 on (— c, c). Now set 

(2.6) T,{VA):=— Ms)iFe)''{s)cosisVA)ds. 



Then each Ts{\/A) is a properly supported smoothing operator and (Te(vA))eg/ 
is an optimal regularization process in the sense of Def. 12.11 From the construction 
of T in terms of the functional calculus of the Laplace operator we conclude the 
following invariance properties of the corresponding embedding lt: 

Corollary 2.5. 

(i) Let f : M ^ M be an isometry. Then for any u G 'D'{M), bT{f*u) — 
f*iT{u). 
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(ii) If "^ is a pseudodifferential operator commuting with A, then Vf commutes 
with It ■ 

For the special case R" with the Euchdean metric the above construction gives 
T,u^fi,^u, with /IJ(e) = -Lj-(0ii?(l))(|^|). 

In particular, for the one-dimensional case n = 1 we obtain fif, = (j)c ■ F{./e)/{2iTe). 
Note that F{./e)/{2TTe) is a standard moUifier, i.e., a Schwartz function with unit 
integral and all higher moments vanishing. Thus the regularization process reduces 
precisely to the usual Colombeau embedding via convolution ( [TOl [SU [20] ) . 

A further consequence of the above construction is that it relates naturally to 
isomorphisms of vector bundles. Thus let : E'l — > £'2 be a vector-bundle isomor- 
phism (covering the identity map on M) and let Di be an order one admissible 
operator on Ei. We choose a Hermitian structure on E2 which makes an isom- 
etry. Then the push-forward operator D2 '■— ipDicf)^^ is isospectral to Di. The 
naturality of functional calculus then gives that for any Schwartz function F with 
F = 1 near the origin the embedding and the bundle map 4> commute, that is 

cf>oF,{Di)=F,{D2)o(f>. 

We note that if 7'i + si ~ r2 + S2 then any Riemannian metric provides an isomor- 
phism of the tensor bundles T^J —> l^^ by 'raising or lowering of indices'. Thus if we 
pick the connection Laplace operators on the tensor bundles and an even Schwartz 
function F as above we obtain a regularization process that is well behaved with 
respect to raising and lowering of indices. 

3. Regularization on globally hyperbolic space-times 

We now return to the situation described at the end of Section[TJ Thus, through- 
out this section, {M,g) will be a smooth space-time, i.e., a connected time-oriented 
Lorentz manifold. We first review the concept of metric splitting for globally hy- 
perbolic space-times. Building on this we construct invariant regularizations of 
distributions on smooth globally hyperbolic space-times. 

3.1. The metric splitting of globally hyperbolic space-times and associ- 
ated Riemann metrics. The original definition, due to J. Leray of global hyper- 
bolicity of a space-time M appeared in ^28j . It requires that the set of causal curves 
connecting two points p, q G M he compact in the space of all rectifiable paths with 
respect to a suitable metric topology ([5J Ch. XII, Sec. 8, 9]). In what follows, we 
will use an equivalent definition, cf. [HI Sec. 6.6], and [51 Ch. XII, Th. 10.2]. 

Thus we call a space-time globally hyperbolic if it satisfies (a) strong causality 
and (b) for any p,q G M the intersection J+ {p) n J~ (q) of the causal future J"*" (p) 
of p with the causal past J~{q) of g is a compact subset of M. Thanks to [3] 
condition (a) may be weakened to causality, i.e., non-existence of closed time- like 
curves. 

As shown by Geroch (cf. [TH]) global hyperbolicity is equivalent to the existence 
of a Cauchy hypersurface and in turn provides a foliation of M by Cauchy hy- 
persurfaces. Further development of these constructions and techniques led to the 
following result due to Bernal- Sanchez (cf. [2|) on the so-called metric splitting of 
a globally hyperbolic space-time {M,g): There exists a Cauchy hypersurface S in 



M and an isometry of {M, g) with the Lorentz manifold (M x S, A) with Lorentz 
metric A given by 

(3.1) X = -Pdt'^ + ht, 

where (3 g C°°(R x S) is positive, (/it)tgE is a smoothly parametrized family of 
Riemannian metrics on S, and t denotes (slightly ambiguously) both the global 
time function (i, x) i— >■ i and its values. In other words, in order to construct a reg- 
ularization (or embedding) for distributions on the globally hyperbolic Lorentzian 
manifold (M, g) we may as well assume that (M, g) = {M. x S, A) . 

The specific structure of the Lorentz metric A in p.ip suggests to associate with 
it the Riemann metric 

(3.2) p:=pdt^ + ht 

on M X 5* and to simply use the regularization and embedding of distributions on M x 
S based on this Riemann metric p. The regularization construction on Riemannian 
manifolds described above requires completeness of the Riemann metric p. It may 
happen that p is not complete, however we may then introduce an appropriate 
conformal factor to obtain a complete Riemannian metric (cf. [30]). Note that the 
latter would amount to introducing the exact same conformal factor for the Lorentz 
metric A and would not change the class in the so-called causal hierarchy of space- 
times according to Minguzzi-Sanchez [23 • Therefore we assume henceforth that p 
is complete. 

Denoting by Ap the Laplacian w.r.t. p, for any s G M and any u G 'D(M) we set 

||u||,:=||(l + Ap)^/2^||i2(M). 

The Sobolev space H^{M) of order s is the completion of 'D{M) with respect to 
this norm. 

Example, (i) Consider 5 = R and A = —dP + if^6^6 as Lorentz metric on R^. It is 
elementary to check that (R^, A) is globally hyperbolic, e.g., by showing that {0} xR 
is a Cauchy hypersurface. The associated Riemann metric p = dt^ + if^6^6 is not 
complete, since the hypersurfaces {t} x M with t y^ are closed and bounded but not 
compact in (M^,p). Multiplying p by the function a G C°°(R2), a{t,x) = l + t^x^, 
yields a conformal metric which is complete (by the Hopf-Rinow theorem, since 
lengths of curves w.r.t. a ■ p are greater than or equal to their Euclidean lengths), 
(ii) Let (S, ho) be a connected Riemannian manifold. For the Robertson- Walker 
space-time R x S" with Lorentz metric of the form A = —dP + f{t)^ho, where 
/ G C°°(R) is positive, we have (cf. [3 Lemma A.5.14]): 

(R X S, A) is globally hyperbolic if and only if (5, hg) is complete. 
If this is the case, then the corresponding Riemannian metric p = dP + f{t)^ho on 
R X S* is complete (see [35], Lemma 7.40). 

3.2. Regularization and embedding via the associated Riemannian struc- 
ture. Let (Tg)^^/ be the regularization (i.e., family of properly supported smooth- 
ing operators) and t: I''(R x S) ^-)- ^(R x S) be the embedding associated with 
the complete Riemann metric p as in (12. 6p . We recall that the embedding thus 
respects the differential algebraic structure of C°°(R x S) and the wave front sets 
of distributions in the strong sense of generalized functions, that is, with respect to 
equality on the level of t/(R x 5*). Furthermore, i is also invariant under isometrics 
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of the Riemannian structure and commutes with the action of the Laplace operator 
Ap corresponding to p on distributions and generaUzed functions, respectively. 

Let u € I?'(R X S). Then according to (I3.ip the difference between the wave 
operator D;^ and the Laplace operator A^ acts on u as 



(3.3) □.. - ^,u . ^=1= 3^'^ 3,u) :^ -28.. 

From this equation and the commutation property of the embedding t with Ap 
we obtain for any u G T>' {W x S) 

(3.4) l{Dxu) - nxL{u) ^ t((Ap - 2e)u) - (Ap - 2e)i(u) = 2(et(u) - l{Qu)) . 

Thus, the precise invariance properties of the embedding with respect to the wave 
operator can be reduced to investigating the corresponding behavior upon inter- 
changing i with Q. 

As a first simple, but useful, observation we point out that t commutes with D\ 
in the distributional sense (i.e., in the sense of association in Q), which follows from 
the fact that Da is a differential operator with smooth symbol. This proves the 
following statement. 

Proposition 3.1. For any u G V {M. x S) we have l{D\u) ~ Da i.{u). 

We will now show that under additional regularity assumptions on the distribu- 
tion u a stronger asymptotic property holds. 

Theorem 3.2. If u E H^{M. x S) is compactly supported then 

Proof. As in (|3.4p we have the basic relation 

Now recalling the Bochner-integral defining the action of T^ on L^-functions we 
have 

QT^u-T^eu^— / (l)c{s)-F{-)(Qcos{sy/-Ap)u- cos{sy/~-Ap)eu\ ds. 

Hence everything boils down to deriving asymptotic estimates for the commutator 
of 9 with the operator S{s) := cos(s^-Ap) (s e K). For any v G L^(M x S) the 
function w G C°°(K,i2(K x 5)), w{s) := S{s)v (s G K), is the mild solution to the 
following Cauchy problem on R x M x S*: 

(*) d^w - Apw = 0, w{0) = V, dsw{{)) = 0. 



In the above integral formula the term cos.{s^J~Ap)Qu corresponds to the solution 
with V = 6m, whereas the term Q cos{syJ~Ap)u is just the application of Q to the 
solution w(s) corresponding to v — u. Applying 6 to (*) we obtain 

dlQw - ApQw = [9, Ap]w =: / 

and 

Qw |.=o= 6(w;(0)) = 9u, dsQw |«=o= Q{dsw{Q)) = 0. 
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Therefore the Duhanicl principle yields 



6cos(s-\/— Ap)u = Qw{s) 

= cos{sy/~Ap)eu+ / (s-r) sinc((s-r)^-Ap)/(r)dr, 
Jo 

where sine: R — >■ R is given by sinc(z) = 2HH£)_ Jj^ summary, we obtain 

QT,u-Teeu^— / Ms)-F{-) / (s-r)sinc((s-r)y^A;)/(r)drds, 

where 

/(r) = [e,Ap]ii;(r) = [6, Ap] cos(ry^A;)u. 

Since [9, Ap] is of third order and maps H^^ into L^ and cos(r-y/— Ap) has 
operator norm 1 on every Sobolev space there exists a constant Ci > such 
that ||/(?')||^2 < C'i||w||/f3 for every r S R. Furthermore, the operator norm of 
(s — r) sinc((s — r)-^/— Ap) is bounded by sup^g^ |(s — r) sinc((s — r)z)\ = \s — r\. 
Combining these upper bounds we estimate 

\\eT,u-T,eu\\^.<^ /'l0^(s)|i|i?(f)l^ds.|l„ll 

Ztt Jt^ eel 

= %\\^\\m I \Me'r)\\Fia)\a' da ■ e' < ^ MnsUch^.WP'h. ■ e\ 

D 

Remark 3.3. Applying the reasoning of the proof of Theorem l3.2l to the operators dt 
or Ma (multiplication by a) instead of O gives the following additional asymptotic 
properties: 

(i) If u e H'^{^ X S) then || T^dtu - dtT.u \\^2 = ©(e^) (e ^ 0). 

(ii) Let a e C°°(R x S) and let u e H^{Rx S) have compact support. Then 

\\T^{au)-aT^u\\^2=0{e^) (e ^ 0) 

We have constructed the regularization operators (Te)eg/ and the embedding t of 
distributions on R x S* based on the Riemannian metric p given by p.2p . However, 
the construction itself does not directly reflect the metric splitting in p.2p or p.ip 
and so far we have not paid special attention to the foliation by the space-like 
Cauchy hypersurfaces {t} x 5 (i G R). However, this foliation becomes essential 
in case of distributions that allow restriction to these hypersurfaces. The latter is 
true in particular for distributional solutions to the wave equation 0\u = /, where 
/ e C°°{R,V'{S)). In these situations we automatically have u £ C°^(R,V'{S)), 
since Char(nA) does not contain any elements of the form (i, x; ±1, 0) G r*(R x 5) 
(cf. [15j, 23.65.5). Thus for every t € R the value u{t) is an element of I?'(S') and the 
metric splitting (|3.ip provides us with a Riemannian metric ht on S, which can be 
used to regularize or embed u{t) according to our general construction. For every 
t g R let {T^^)e£i denote the regularization obtained from the Riemann metric ht 
on S and let tht denote the corresponding embedding V^S) '-> G{S). The following 
statements compare these with the global constructions on R x S*. 

Theorem 3.4. (i) The embedding l: !)'(& x S) ^^ G{M. x S) respects the metric 
splitting (|3.ip in the weak sense: If u £ C°°(R, I''(S')) then we have 

Vie R: iiu){t) « iht{u{t))- 
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(a) If u £ C°° {M., H'^ (S)) , then we have for every compact subset Z <ZR 

sup II {T,u){t) - T^^u{t) ll^.(^) = 0{e^) {e ^ 0). 
tez 

Proof, (ii) Pick x e 2?(M) such that x = 1 near Z and let u{t) := x{t)u{t). It then 
foUows from [13], Prop. 3.7 that {Ti^u — T^u)e is negligible on supp(u — uY 5 Z . 
Thus we may without loss of generality assume that the support of u is bounded 
in i, so u G L^(K x S). This allows us to employ the integral formulae defining Tg 
and Tg''* and obtain 

{T,u){t)-T^*u{t) = 

^ [ (/)e(s)-F(-)( (cos(sv/^A;)w)(t) - cos{sy/~AhJu{t) ) ds, 

^'^ JR £ £ Vv ^ ' V ^ '/ 

— -.w{s,t) —■.v{s,t) 

where v,w E C°°(IR^,i?^(S')) solve the following Cauchy problems, respectively: 

{d^,~Ap)w^O (mV'iRxRxS)), w \s=o = u, dsw\s=o = Q, 

{dl - Ahjv{., t) = (in I?'(R x S)), v{0, t) = u{t), dsv{0, t) = 0. 

Therefore we have at arbitrary, but fixed t G M 

dlw{.,t) - AhM;t) = Apu;(.,t) - AhM;t) =■ ft, 

where /( £ C°°(M., L'^{S)), since u (hence w) has values in H^{S). We observe that 
tpt '■= w{.,t) — v{.,t) satisfies the Cauchy problem 

(af - Aft JV^t = /* (in I?'(R X 5)), M0) = 0, 9.Vt = 0, 

which implies 

'^t{s) = / (s-r) sinc((s - r)^J-Ah,)ft{r)dr. 
Jo 



Hence we obtain 



S2 



||V't(s)llL2(5) < I / (s-r)||/t(r)||^2(5)dr| < — • sup |l/t(cr)||L2(s) 
JO ^ l<'"l<kl 

and therefore arrive at 



\\{T,u){t)~T,^'uit)\\^,^,^<— \^^{s)\-\F{-)\\\M^)\\^,^s) 

^^ JR £ £ 

<^ r^ \UeT)\\F{T)\T^ sup \\ft{cT)\\I^.^s)dT 
47r J-2c/e W\<e\T\ ' 



£ 



2 



<^Uc\\l^ sup ||/*(a)||^.(5)||^"||^, =0(e^) (e ^ 0) 

uniformly when t varies in a compact set. 

(i) Since association is checked by action on test functions in 'D{S) we may reduce 
to the case that supp(w(t)) is compact. Then u{t) £ H^{S) for some I £ Z and since 
(1 — A/ij)^'/^ commutes weakly with i and strongly with ih^ the assertion follows 
from (ii). D 
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